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Abstract 

Flavour inclusive, udsc and b fragmentation functions in unbiased jets, and flavour 
inclusive, udsc, b and gluon fragmentation functions in biased jets are measured 
in e + e~ annihilations from data collected at centre-of-mass energies of 91.2, and 
183-209 GeV with the OPAL detector at LEP. The unbiased jets are defined by 
hemispheres of inclusive hadronic events, while the biased jet measurements are 
based on three-jet events selected with jet algorithms. Several methods are employed 
to extract the fragmentation functions over a wide range of scales. Possible biases 
are studied in the results obtained. The fragmentation functions are compared to 
results from lower energy e + e~ experiments and with earlier LEP measurements 
and are found to be consistent. Scaling violations are observed and are found to be 
stronger for the fragmentation functions of gluon jets than for those of quarks. The 
measured fragmentation functions are compared to three recent theoretical next- 
to-leading order calculations and to the predictions of three Monte Carlo event 
generators. While the Monte Carlo models are in good agreement with the data, 
the theoretical predictions fail to describe the full set of results, in particular the b 
and gluon jet measurements. 
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1 Introduction 

Hadron production in high energy collisions can be described by parton showers (successive 
gluon emissions and splittings), followed by the formation of hadrons which cannot be described 
perturbatively. Gluon emission, the dominant process in parton showers, is proportional to the 
colour factor associated with the coupling of the emitted gluon to the emitter. These colour 
factors are Ca = 3 when the emitter is a gluon and Cf = 4/3 when it is a quark. Consequently, 
the multiplicity of soft gluons from a gluon source is (asymptotically) 9/4 times higher than 
from a quark source [1]. The inequality between Ca and Cf plays a key role in the explanation 
of the observed differences between quark and gluon jets: compared to quark jets, gluon jets 
are observed to have larger widths [2], higher multiplicities [2,3], softer fragmentation func- 
tions [2,4,5], and stronger scaling violations of the fragmentation functions [5]. 

The fragmentation function, D^(x,Q 2 ), is defined as the probability that parton a, which 
is produced at short distance, of order 1/Q, fragments into hadron, h, carrying fraction x of 
the momentum of a. In this study, the momentum fraction is defined as Xe = Eh/Ej et , where 
Eh is the energy of the hadron h and Ej et is the energy of the jet to which it is assigned. The 
relative softness of the gluon jet fragmentation function is explained in the low xe region by 
the higher multiplicity of soft gluons radiated, and in the high xe region by the fact that the 
gluon cannot be present as a valence parton inside a produced hadron (first a splitting g — > qq 
has to occur). The stronger scaling violation is due to the fact that the scale dependence of 
the gluon jet fragmentation function is dominated by the splitting function P g ^ gg ~ Ca, while 
that of the quark jet is dominated by the splitting function -Pq-+qg ~ Cf- 

Jets in e + e~ annihilations are commonly defined using a jet finding algorithm, which is a 
mathematical prescription for dividing an event into parts associated with individual quarks 
and gluons. For example, quark and gluon jets are often defined by applying a jet finder to 
select three-jet qqg events. Some of the most common algorithms are the Durham [6] and 
cone [7] jet finders. Different jet finders result in different assignments of particles to jets: thus 
jets defined using a jet finding algorithm are called biased. In contrast, quark and gluon jets 
used in theoretical calculations are usually based on inclusive samples of back-to-back qq and 
gg final states rather than three-jet events. A hemisphere of a qq event is defined as a quark 
jet and similarly, a gluon jet is defined by a hemisphere in a gg final state. The hemisphere 
definition yields a so-called unbiased jet because the jet properties do not depend on the choice 
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of a jet finder. Measurements of unbiased quark jets have been performed at many scales since 
such jets correspond to hemispheres of inclusive e + e~ — > hadrons events [8-10]. Direct mea- 
surements of unbiased gluon jets are so far available only from the CLEO [11] and OPAL [4,12] 
experiments, however. At CLEO, jets originating from radiative T decays have energies of only 
about 5 GeV, which limits their usefulness for jet studies. In [4, 12], unbiased gluon jets were 
selected using rare events of the type e + e~ — > qqgmcb i n which the object g inc i, taken to be the 
gluon jet, is defined by all particles observed in the hemisphere opposite to that containing two 
b-tagged quark jets. Due to the low probability of such a topology, this method of obtaining 
unbiased gluon jets is only viable for very high statistics data samples. Recently, the OPAL 
experiment has measured properties of unbiased gluon jets indirectly. In [13], recent theoretical 
expressions to account for biases from event selection were used to extract gluon jet properties 
over a range of jet energies from about 11 to 30 GeV. In [14], the first experimental results 
based on the so-called jet boost algorithm, a technique to select unbiased gluon jets in e + e~ 
annihilations, were presented for jet energies from 5 to 18 GeV. 

Scaling violations of quark and gluon jet fragmentation functions from three-jet events pro- 
duced in e + e~ collisions at a center-of-mass system (c.m.s.) energy of \fs = 91.2 GeV, based on 
the kx jet algorithms Durham [6] and Cambridge [15], were reported in [5]. These scaling vio- 
lations were found to be consistent with the expectations from the Dokshitzer-Gribov-Lipatov- 
Altarelli-Parisi (DGLAP) evolution equations [16]. In our study, we present measurements of 
quark and gluon jet fragmentation functions at y/s = 91.2 GeV and y/s = 183-209 GeV. The 
data were collected with the OPAL detector at the LEP e + e~ collider at CERN. We measured 
seven types of fragmentation functions: the udsc, b, gluon and flavour inclusive fragmenta- 
tion functions in biased jets, and the udsc, b, and flavour inclusive fragmentation functions 
in unbiased jets. While the two types of flavour inclusive fragmentation functions have been 
measured many times, data on the other types of fragmentation functions are still rather scarce. 

The paper is organised as follows. In Section 2, a brief description of the OPAL detector is 
given. The samples of data and simulated events used in the analysis are described in Section 3. 
In Section 4, the event and jet selections are discussed. The analysis procedure, including the 
methods used to evaluate systematic uncertainties, is presented in Section 5. Section 6 deals 
with a Monte Carlo (MC) study of the biases introduced by our jet finding procedure. Next-to- 
leading order (NLO) calculations [17-19] for fragmentation functions are described in Section 7. 
In Section 8, we present a comparison of our data to other measurements, to MC predictions, 
and to the NLO calculations. A summary and conclusions are given in Section 9. 

2 The OPAL detector 

The OPAL detector is described in detail elsewhere [20]. The tracking system consists of a 
silicon microvertex detector, an inner vertex chamber, a large volume jet chamber and spe- 
cialized chambers at the outer radius of the jet chamber which improve the measurements in 
the ^-direction 1 . The tracking system covers the region |cos6>| < 0.98 and is enclosed by a 
solenoidal magnet with an axial field of 0.435 T. Electromagnetic energy is measured by a 
lead-glass calorimeter located outside the magnet coil, which covers | cos#| < 0.98. 

1 OPAL uses a right-handed coordinate system defined with positive z along the electron beam direction and 
with positive x pointing towards the centre of the LEP ring. The polar angle 9 is defined relative to the +z 
axis and the azimuthal angle relative to the +x axis. 
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3 Data and Monte Carlo samples 



The present analysis is based on two data samples which will be referred to as the LEP1 and 
LEP2 samples. The LEP1 data sample contains hadronic Z decay events collected with the 
OPAL detector between 1993 and 1995 at c.m.s. energies within 250 MeV of the Z peak. The 
LEP2 data sample contains hadronic events collected with the OPAL detector in the period 
1997-2000 at c.m.s. energies in the range 183-209 GeV. All the data were taken with full 
readout of the r-0 and z coordinates of the silicon microvertex detector which is essential for 
precise measurements of primary and secondary vertices. The total integrated luminosity in 
the LEP1 data is 130 pb _1 , while the LEP2 data sample corresponds to a luminosity of 690 pb _1 . 

In this study, we work with three types of MC event samples. The detector level samples 
include full simulation of the detector response [21], the initial-state photon radiation (ISR) 
and background processes, and contain only those events which pass the same selection cuts as 
applied to the data. The hadron level samples do not include ISR or detector simulation and 
allow all particles with lifetimes shorter than 3 xlO" 10 s to decay. The parton level samples 
are formed by final-state partons, i.e. quarks and gluons present at the end of the perturbative 
shower, and do not include ISR. 

Signal MC events for the LEP1 data, of the form e + e~ — > Z — > qq(g), are generated using 
the JETSET 7.4 [22] and HERWIG 6.2 [23] programs with the parameter settings tuned on 
LEP1 OPAL data described in [24] and [25], respectively. For LEP2 data, the signal e + e~ — > 
Z*/ 7 * _> qq( g ) events are simulated using PYTHIA 6.125 [22,26] and HERWIG 6.2. For events 
of this type, PYTHIA is the same as JETSET except for the inclusion of ISR processes. The 
same parameter settings are used for the LEP2 PYTHIA and HERWIG samples as are used for 
the LEP1 JETSET and HERWIG samples, respectively. In the JETSET and PYTHIA event 
generators, the string fragmentation model is implemented, while HERWIG uses the cluster 
fragmentation model. The initial- and final-state photon radiation for the LEP2 MC samples 
are performed by interfacing the KK2F program [27] to the main generator programs. In addi- 
tion to PYTHIA and HERWIG we also use the ARIADNE 4.08 [28] event generator to compare 
with the final results. For hadronization, the generator is interfaced to the JETSET 7.4 pro- 
gram. The parameter settings used for ARIADNE are documented in [4,29]. 

To estimate the background in the LEP2 data, we generate events of the type e + e~— > 
4 fermions. These events, in particular those with four quarks in the final state, constitute 
the major background in this analysis. The 4-fermion events are generated using the GRC4F 
2.1 [30] MC event program. The final states are produced via s-channel or t-channel diagrams 
and include W + W~ and ZZ events. This generator is interfaced to PYTHIA using the same 
parameter set for the fragmentation and decays as used for the signal events. 

The signal as well as the background MC event samples for the LEP2 period are generated 
at c.m.s. energies of 183, 189, 192, 196, 200, 202, 204, 205, 206, 207 and 208 GeV reflecting the 
energy distribution in the collected data samples. 
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4 Event and jet selection 



4.1 Selection of hadronic Z and Z*/7* events 

The procedures for identifying hadronic events are discussed in [31]. The selection of the inclu- 
sive hadronic event sample in the LEP1 data is based on tracks and electromagnetic clusters. 
Tracks are required to have at least 40 measured points (of 159 possible) in the jet chamber, to 
have a momentum greater than 0.15 GeV/c, to lie in the region | cos#| < 0.94, to have a distance 
of the point of closest approach to the collision point in the r-<p plane, d < 5 cm, and along 
the z axis, z$ < 25 cm. Clusters are required to be spread over at least two lead glass blocks 
and to have an energy greater than 0.10 GeV if they are in the barrel section of the detector 
(| cos 6*| < 0.82) or greater than 0.20 GeV if they are in the endcap section (0.82 < | cos6>| < 0.98). 
A matching algorithm is employed to reduce double counting of energy in cases where tracks 
point towards electromagnetic clusters. Specifically, the expected calorimeter energy of the 
associated tracks is subtracted from the cluster energy. If the energy of a cluster is smaller 
than that expected for the associated tracks, the cluster is not used. Each accepted track and 
cluster is considered to be a particle. Tracks are assigned the pion mass. Clusters are assigned 
zero mass since they originate mostly from photons. 

To eliminate residual backgrounds, the number of accepted tracks in each event is required 
to be at least five. To reject events in which a significant number of particles is lost near the 
beam line direction, the thrust axis of the event, calculated using the particles, is required to 
satisfy | cos(0 t hrust)| < 0.90, where thrust is the angle between the thrust and beam axes. The 
two-photon background (events of the type 77 — * qq) is reduced by imposing the conditions 
-Evis/v^ > 0.1 and |pbai| < 0.6, where E v[s is the total visible energy (i.e. the sum of de- 
tected particle energies) and pbai is the momentum sum in the z direction, normalized by E v - ls . 
The residual background in the LEP1 data sample from all sources is estimated to be less 
than 1% [31]. The number of inclusive hadronic events is 2 387227 (see the first row in Table 
1), with the selection efficiency estimated to be 96%. 

At c.m.s. energies above the Z resonance, several new sources of background exist. To se- 
lect hadronic events in the LEP2 data, the same procedure as described for the LEP1 data is 
used and in addition, we apply the procedure described in [32-35] to reduce the background as 
summarized below. 

The majority of hadronic events at LEP2 are radiative events in which initial-state radiation 
reduces the original c.m.s. energy of the hadronic system. To reject such ISR events, we deter- 
mine the effective c.m.s. energy of the hadronic system, \fs i , following the procedure described 
in [35] which takes possible multiple photon radiation into account. We require y/s — \fs i < 10 
(20) GeV to select inclusive hadronic events for the hemisphere (three-jet) analysis described 
below. We refer to this procedure as the "invariant mass" selection. For systematic studies, we 
apply an alternative method based on combining cuts on the visible energy and missing mo- 
mentum of the event and on the energy of an isolated photon candidate [32]. This procedure 
is referred to as the "energy balance" selection. Simulated hadronic Z*/7* events are defined 
to be radiative if s' tiue < y/s — 1 GeV, where s[ rue is the true effective c.m.s. energy. The 
efficiency for selecting LEP2 non-radiative hadronic events is 73%. 

The production of W + W~ and ZZ pairs with hadronic or semi-leptonic decays (4-fermion 
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final states) is an additional source of background. This background is reduced by applying 
a method described in [35]: first each event is forced into a four-jet configuration using the 
Durham jet finder. In the LEP2 samples, the 4-momenta of all measured particles are boosted 
into the rest frame of the hadronic system with the effective c.m.s. energy, v^s 7 , and are then 
used to find jets. Then an event weight Wqqb is defined based on calculated QCD matrix 
elements for the process e + e~ — > qqqq or qqgg, with the four parton final state corresponding 
to the obtained four-jet kinematics [36]. The QCD matrix elements are calculated using the 
EVENT2 program [37]. A good separation between the Z*/7* and W + W~ or ZZ pair events is 
achieved by requiring Wqcd > —0.5. 

The remaining background from e + e _ — > t + t~ and two-photon events is estimated to be 
about 0.2% [35] and is neglected. The remaining 4-fermion background is subtracted from the 
data bin-by-bin. The number of the inclusive hadronic events in the LEP2 data sample for the 
hemisphere (three-jet) analysis is 10 866 (12 653) with 11% (14%) 4-fermion background (see 
the first row in Table 1). 

4.2 Jet selection 

As explained in the introduction, we employ two definitions of jets. In the inclusive hadronic 
event samples we use the unbiased jet definition where the jets are defined by particles in 
hemispheres of the qq system. In the three-jet samples, we apply a jet algorithm and thus 
work with biased jets. Three jet algorithms are used: the Durham [6], Cambridge [15] and 
cone [7] algorithms. Relatively large differences in the techniques used by the kr and cone jet 
finders ensure that the jet finder dependence of the results is estimated conservatively. The 
jet algorithm is forced to resolve three jets in each event. The jet energies and momenta are 
then recalculated by imposing overall energy-momentum conservation with planar massless 
kinematics, using the jet directions found by the jet algorithm. The jet energies are given by 
the relation: 

Ei = V ~ • °' M (1) 

sin $ij + sin 9 j:k + sin 6 kyi 

where 6^ is the angle between jets i and j and k corresponds to the remaining jet. We note 
that for the LEP2 detector level jets, the effective c.m.s. energy, V^, is used in the above 
formula. Cuts, given in Table 2, are chosen to ensure that the jets are well contained within 
the sensitive part of the detector, well separated from each other and that the event is planar. 
The numbers of LEP1 and LEP2 events passing these selection criteria are shown in the second 
row of Table 1. The efficiency for selecting non-radiative three-jet LEP2 events is 68%. 

All three jet algorithms yield very similar jet angular and energy resolutions, with the 
Durham algorithm being slightly better than the other two. Therefore, the Durham algorithm 
is used as the reference, with the cone and Cambridge jet finders used for systematic studies. 
The jet energy resolution, defined as (£g tector - Ef^™) / Ef^ on , is found to range from 2% for 
the most energetic jet to 11% for the least energetic jet. The distribution of the angles between 
the detector and parton jet axes is found to have an RMS of 0.05 radians for the most energetic 
jet and 0.16 radians for the least energetic jet. See Section 5.2.1 for an explanation of how the 
detector and parton level jets are associated with each other. 
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5 Analysis procedure 



In the following, we describe the method we use to determine the quark and gluon jet fragmen- 
tation functions. The measured fragmentation function is defined here as the total number of 
charged particles, N p , in bins of x E and scale Q normalized to the number of jets, Nj et (Q), in 
the bin of Q: 

J- dN p (x E , Q) 

N ieb (Q) dx E 1 ' 

where x E is defined in the Introduction. 

5.1 Jet scale Qj et 

To measure the scale dependence, it is necessary to specify a scale relevant to the process under 
study. For inclusive hadronic events, the scale is y/s. For jets in three-jet events, neither y/s 
nor £j e t is considered to be an appropriate choice of the scale [38]. QCD coherence suggests [39] 
that the event topology (i.e. the positions of the partons with respect to each other) should 
also be taken into account. In studies of quark and gluon jet characteristics [5,38,40,41] the 
transverse momentum-like scale Qjet, of a jet with energy £j et has been used: 

•d 

Qjet = Ejetsin-, (3) 

where d is the angle between this jet and the closest other jet. This scale roughly expresses 
a maximum allowed transverse momentum (or virtuality) of gluons radiated in the showering 
process with respect to the initial parton, whilst still being associated with the same jet. This 
definition of scale is adopted for the present analysis. The jet energy and Qjet spectra are shown 
in Figs. 1 and 2 for the three jets found by the Durham jet algorithm and ordered in energy, 
with jet 1 being the most energetic and jet 3 the least energetic jet. The data are seen to be 
well described by the JETSET and HERWIG models. A similar description is also seen for the 
cone and Cambridge jet finders (not shown). 

5.2 Quark and gluon jet identification 

There are several ways to identify quark and gluon jets. In this analysis, three methods are used: 
the b-tag and the energy-ordering methods to identify quark and gluon jets in biased three-jet 
events, and the hemisphere method to identify unbiased quark jets in inclusive hadronic events. 
In addition, b tagging is used to separate udsc and b quark jets from each other, both for 
the biased and unbiased jet samples. In contrast to the b-tag method, the energy-ordering 
method only allows flavour inclusive quark jets to be distinguished from gluon jets. Note that 
the flavour composition of the primary quarks in e + e~ — > qq is predicted by electroweak theory 
to vary with c.m.s. energy. Therefore, to perform a meaningful comparison of the biased jet 
data taken at y/s = 91.2 GeV with the unbiased jet data measured at several c.m.s. energies, a 
special correction is applied in the construction of the flavour inclusive fragmentation function 
from biased jets (see Section 5.3). 

5.2.1 b-tag method in three-jet events 

In the three-jet sample, the b-tagging technique is used to obtain samples enriched in udsc, b or 
gluon jets. The analysis utilizes an inclusive single jet tag method based on a neural network, 
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as described in [42]. Any or all of the three jets may be used to extract the fragmentation 
functions. Note that with our selection of three-jet events, the highest energy parton jet is 
predicted to be the gluon jet in 4.8% of the events. 

In the data and MC, three samples of jets are selected, each with different fractions orig- 
inating from udsc-quarks, b-quarks or gluons. We first look for jets with secondary vertices 
found in cones of radius R = 0.65 radians from the jet axes. A jet is considered to be a b-tag 
jet if it contains a secondary vertex with neural network output value, VNN, greater than 0.8 
for LEP1 events or 0.65 for LEP2 events. A jet with no secondary vertex, or a vertex with 
VNN<0.5 is considered to be an "anti-tag" jet. The b-tag and gluon jet samples are taken 
from events with one or two b-tag jets and at least one anti-tag jet. If one or two b-tag jets 
and one anti-tag jet are found, the b-tag jets enter the b-tag jet sample and the anti-tag jet 
enters the gluon jet sample. If one b-tag jet and two anti-tag jets are found, the b-tag jet enters 
the b-tag jet sample, and the lower energy other jet is included in the gluon jet sample. The 
udsc jet sample is formed by all three jets in events with no b-tag jet or with b-tag jets but no 
anti-tag jet (the contribution from the latter events is negligible in practice). Note that with 
this definition, the gluon jet is explicitly included in the udsc jet sample. The correction proce- 
dure to obtain a pure udsc jet sample with the gluon jet component removed is described below. 

The purities of the different jet samples are evaluated by examining Monte Carlo events 
at the parton, hadron and detector levels. First, parton level jets are examined to determine 
whether they originate from a quark or a gluon. This determination is performed in two ways: 

• Flavour assignment: It is assumed that the highest momentum quark and antiquark 
with the correct flavour for the event are the primary quark and antiquark. In events in 
which different parton level jets contain the primary quark and antiquark, the remaining 
jet is assumed to arise from a gluon. 

• Non-flavour assignment: A parton jet is identified as a quark (antiquark) jet if it 
contains an arbitrary number of qq pairs and gluons plus one unpaired quark (antiquark). 
If such two parton jets are found, the gluon jet is defined as that containing only qq pairs 
(if any) and gluons. 

A small fraction of events showing an ambiguous assignment of the primary qq pair and gluon 
to three parton level jets is excluded from the event samples. It amounts to 1.3% for the flavour 
and 2.5% for the non-flavour assignment. To obtain the final results, the former method is used. 

Detector and parton level jets are assigned to the hadron jet to which they are nearest in 
angle. For events in which more than one parton or detector level jets are assigned to the same 
hadron level jet (about 9% of the events), the closest jet is chosen, while the more distant jet 
is assigned to the remaining hadron jet. The above procedure is referred to as the "matching" 
procedure, and the hadron level jets associated with the parton level quark and anti-quark jets 
are defined to be pure quark jets, while the remaining jet is a pure gluon jet. 

The purity and the efficiency of the LEP1 and LEP2 b-tag jet samples as a function of the 
VNN variable are shown in Fig. 3. The purity of the b-tag jet sample at the point VNN=X is 
defined as the fraction of pure b jets in the sample of b-tag jets with VNN>X. The efficiency 
of the b-tag jet sample at the point VNN=X is defined as the fraction of the b-tag jets with 
VNN>X in the sample of all pure b jets. For VNN>0.8 applied in the LEP1 samples, the purity 
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of the b-tag jet sample is 90% and the efficiency 23%. The corresponding gluon jet purity and 
efficiency are 84% and 40%, respectively. The LEP2 samples are treated analogously to the 
LEP1 samples, except that we require VNN>0.65 because of low event statistics. The b (gluon) 
jet tagging efficiency is 27% (45%) and the purity 60% (80%). 

To obtain a distribution of a variable D (e.g. the fragmentation function) of pure udsc (b, 
gluon) jets, -Df^g), one has to solve the following equation 

uncor pure 

/ Pn flb Pi g \ / A \ 
x E ,Q) = PbiPbbPbg (Q) A \(xe,Q) (4) 
V Pgj P gb P sg J \D g J 

where A^g) s t an ds for a distribution of the variable D obtained from the sample of detector 
level udsc (b-tag, gluon) jets. The purity P^ denotes the probability that a jet from the jet 
sample i comes from a parton j. The indices i, j run over symbols l,b and g which stand for the 
u,d,s,c ( "light" ) -quark, b-quark and gluon. 

In Fig. 4 the LEP1 and LEP2 purity matrices as functions of Qj e t are shown as obtained 
using the Durham jet algorithm. The numbers of selected udsc, b-tag and gluon jets are shown 
in Table 1. The larger number of b-tag jets compared to gluon jets is due to the inclusive single 
jet tag method which allows up to two b-tag jets per event. 

5.2.2 Energy-ordering method 

This method is based on the QCD prediction that in a three-jet event the lowest energy jet 
has the highest probability to arise from a gluon. In this method only jets 2 and 3 are used, 
which form the quark and gluon jet samples, respectively. There are two ways of estimating 
the purities: either via the matching which employs the inter-jet angles as described in the 
b-tag method, or using matrix elements. It has been shown [43] that, for leading order QCD 
matrix elements, the probability for a given jet % among the jets {i, j, k} to be a gluon jet can 
be expressed as a function of the jet energies: 

1§0C (l-xj)(l-x k ) [b) 
where Xi = 2Ej etji /y/s. The corresponding probability for the jet to be a quark jet is 

= 1 - Pig, (6) 

normalised such that P lq + P2 q + Ps q = 2. Thus, in this way, the purities can be obtained based 
on the kinematics of the data, without recourse to MC information. The scale dependence of 
the quark purities of jets 1 and 2, and the gluon purity of the jet 3, are shown in Fig. 5. Good 
agreement is obtained between the data and MC for the matrix element method. The MC 
results based on matching are seen to agree well with the results based on the matrix elements. 
For consistency reasons, the purities based on the matching are used to obtain the final results. 

An unfolding to the level of pure quark and gluon jets is carried out by solving the following 
equation: 

uncor pure 
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where D^" r is the detector level distribution of a variable D in the sample of jets 2 (3) and 
^q(g) C corres P on ds to pure quark (gluon) jets. The energy-ordering method can only be applied 
in the Qj e t region where the samples of jets 2 and 3 overlap (6 < Qj e t < 27 GeV for the LEP1 
and 10 < Q jet < 60 GeV for the LEP2 sample). 

5.2.3 Hemisphere method 

In the inclusive hadronic event sample we again use b-tagging to obtain samples enriched in b 
and udsc jets. In the LEP1 sample, a b-tag event is defined by requiring two secondary vertices 
with VNN>0.8, while in the LEP2 sample — due to limited statistics — only one secondary vertex 
with VNN>0.8 is required. All remaining events form the udsc event sample. Events with no 
requirement on the presence of a secondary vertex form the inclusive hadronic event sample. 
Each event contains two unbiased jets (hemispheres) of the same energy, y/s/2. The jets are 
unfolded to the level of pure udsc and b jets using an analogous procedure to that described 
in Section 5.2.2 for the energy-ordering method. In Eq. (7), we replace the indices 2 and q 
by the index 1, and the indices 3 and g by the index b. The purity P b i (Pi) then denotes the 
probability that a jet from the b-tag (udsc) jet sample comes from an u,d,s or c-quark. In 
the LEP1 MC sample, Pi =79% and P bb =99.7% which means that we work with a very pure 
b-tag jet sample. The corresponding purities for the LEP2 MC sample are 89% and 75%. The 
numbers of unbiased udsc and b-tag jets passing the selection cuts together with background 
estimates (for the LEP2 data) are summarized in Table 1. The higher efficiency of selecting 
non-ISR events and the lower background compared to those for the biased jets is due to the 
tightened cut on the c.m.s. energy described in Section 4.1. 

5.3 Construction of flavour inclusive fragmentation function from 
biased jets 

To construct the flavour inclusive fragmentation function from the LEP1 biased jets, the samples 
of udsc, b-tag and gluon jets from the b-tag method are used. The quark jet sample is formed 
by a sum of the udsc and b-tag jet samples. The unfolding to the level of pure quark and 
gluon jets can then proceed by use of Eq. (7) where the sample of jets 2 is replaced by the 
quark jet sample and the sample of jets 3 by the gluon jet sample. To take into account the ^/s 
dependence of the flavour composition of the primary qq pair, the sample of pure quark jets 
is constructed as a sum of samples of pure udsc and b jets, weighted by factors of r u d sc ((Qjet)) 
and r b ((Qj Ct )), respectively. The r b ((Qj e t)) factor is calculated using the hadron level MC, as 
the ratio of the bb production rate for a given Qj et bin with a mean value of (Qjet) m three-jet 
events generated at y/s = 91.2 GeV and the bb production rate in inclusive hadronic events 
generated at ^/s = (Qjet)- The factor r u d sc ((Qjet)) is determined in an analogous fashion. The 
corrections based on r b and r udsc are smaller than 15% and bring the biased jet data closer to 
the published unbiased jet data. 

5.4 Correction procedure 

The remaining 4-fermion background in the LEP2 data is estimated for each observable by 
MC simulation and subtracted on a bin-by-bin basis from the data distributions, as already 
mentioned in Section 4.1. Then the data and MC distributions at the detector level are corrected 
to the level of pure quarks and gluons by solving either Eq. (4) or Eq. (7). As a last step, we 
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correct the data for the effects of limited detector acceptance and resolution as well as for 
the presence of remaining radiative events. The data are multiplied, bin-by-bin, by correction 
factors calculated as ratios of distributions at the hadron level to those at the detector level. 
For the hadron level biased jets, the same jet selection criteria as described in Section 4.2 
are applied except that the jets are not required to satisfy |cos#j C t| < 0.90. The quark and 
gluon jets at the hadron level are identified with MC information using the matching technique 
described in Section 5.2.1. The correction factors from JETSET/PYTHIA used to correct the 
data do not exceed 20%. The correction factors from HERWIG used to estimate the model 
dependence of the results are similar. A bin-by-bin correction procedure is suitable for the 
measured distributions as the detector and ISR effects do not cause significant migration (and 
therefore correlation) between bins. Typical bin purities for the Qj et binning chosen were found 
to be 75%, the lowest value was 65%. 

5.5 Systematic uncertainties 

The systematic uncertainties of the measurements are assessed by repeating the analysis with 
the following variations to the standard analysis. 

1. The systematics on the modelling of the Z and Z*/7* events used to correct the data 
for ISR, detector effects and quark and gluon jet misidentification is estimated by using 
HERWIG instead of JETSET/PYTHIA. In the bulk of the measured data, the maximum 
differences for all types of fragmentation functions do not exceed 6%. In the last xe bin of 
both types of flavour inclusive fragmentation functions (0.8 < x E < 1.0), the two models 
deviate from each other by as much as 50-60%. 

2. To assess any inadequacies in the simulation of the response of the detector in the endcap 
regions, the analysis was restricted to the barrel region of the detector, requiring the tracks 
and electromagnetic clusters to lie within the range | cos # P articic | < 0.70. The maximum 
differences reach 10% for biased jets (for large x E ) and 2% for unbiased jets. 

3. Potential sensitivity of the results to details of the track selection is assessed by repeating 
the analysis with modified track selection criteria: the maximum allowed distance of the 
point of closest approach of a track to the collision point in the r — (f) plane, do, is changed 
from 5 to 2 cm, the maximal distance in the z direction, z , from 25 to 10 cm and the 
minimal number of hits from 40 to 80. The quadratic sum over the deviations from the 
standard result, obtained from each of these variations, is included to the total systematic 
uncertainty. In most of the bins, the changes are below 1%. Larger changes are observed 
for high xe, where they are within 7% for both, the biased and unbiased jets. 

4. The jet algorithm dependence of the biased jet results is estimated by repeating the 
analysis using Cambridge and cone jet algorithms. The largest of the two deviations from 
the standard result (the cone algorithm in most of the bins) is taken as the systematic 
uncertainty. All differences are within 10% for all types of fragmentation functions, except 
at low Qjet and xe (4 < Qjet < 9 GeV with 0.02 < x E < 0.04) where the results of the cone 
algorithm are about 20%, 24%, 31% and 36% below the results of the Durham algorithm, 
for the flavour inclusive, udsc, b and gluon jet fragmentation function, respectively. The 
differences between the results for individual jet algorithms diminish with increasing jet 
energy. 
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5. The jet selection criteria were varied. The minimum particle multiplicity per jet is changed 
from 2 to 4; the minimum corrected jet energy is changed from 5 GeV to 3 and 7 GeV; 
the minimum inter-jet angle is changed from 30° to 25° and 35° and the minimum sum of 
inter-jet angles is changed from 358° to 356° and 359°. The largest deviation with respect 
to the standard result is taken as the systematic uncertainty. The differences are below 
2% in all cases, except for large x E with small Q jet where they reach 6%. 

6. The dependence of the results on the neural network output value is estimated by varying 
the cut on VNN from 0.50 to 0.95. The maximum of the deviations with respect to the 
standard result is taken as the systematic uncertainty. Typical deviations are 2% for 
unbiased jets and the LEP1 biased jets, while they are 5% for the LEP2 biased jets. The 
largest deviation is 11% for the unbiased jets and 20% for the biased jets (both observed 
for large xe). 

7. The b-tagging efficiency is determined using MC events. The systematic uncertainty in 
this efficiency was estimated to be about 5% for VNN>0.50 in LEP2 data [44]. The effect 
of this uncertainty is assessed by changing the VNN thresholds in the MC samples such 
that the b-tagging efficiency increases or decreases by 10%, while leaving the thresholds 
in the data unchanged. The largest deviation with respect to the standard result is taken 
as the systematic uncertainty. In most of the bins, the differences are below 1%. In the 
high xe region, they reach 4% for unbiased jets and are typically within 8% for biased 
jets. 

8. The uncertainty in the estimates of purities for the b-tag method is accounted for by using 
the non-flavour assignment instead of the flavour assignment of the outgoing primary qq 
pair and gluon to three parton jets. Non-negligible differences in the purities are seen only 
in those Qj et regions where the purities are small. This results in negligible effects on the 
final results: they are below 1% everywhere. In case of the energy-ordering method, the 
procedure based on the matrix elements is used instead of the matching. The differences 
for the gluon jet fragmentation functions are below 1% everywhere. 

9. Uncertainties arising from the selection of non-radiative LEP2 events are estimated by 
using the "energy balance" procedure instead of the "invariant mass" procedure. The 
differences are below 5% for both the biased and unbiased jets. 

10. Systematic uncertainties associated with the subtraction of the 4-fermion background 
events in the LEP2 samples are estimated by varying the cut on W^qcd from -0.5 to 0.0 
and -0.8. The maximum of the deviations with respect to the standard result is taken 
as the systematic uncertainty. The differences are below 4% for both the biased and 
unbiased jets. In addition, we varied the predicted background to be subtracted by ±5%, 
slightly more than its measured uncertainty at y/s = 189 GeV of 4% [45]. The differences 
are below 1% everywhere. 

The results for the udsc jets are found to be less sensitive to the above variations than the 
results for b and gluon jets. The largest changes in the numbers of selected b and gluon jets 
relative to those shown in Table 1 are given by variation 6. For the LEP1 sample, the number of 
b-tag (gluon) jets grows by 55% (44%) for VNN=0.5 and drops by 40% (37%) for VNN=0.95. 
Variation 6 also gives rise to the most significant change in the purities of the b-tag and gluon 
jet samples. The b (gluon) purity decreases by 17% (5%) for VNN=0.5, while it increases by 
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7% (2%) for VNN=0.95 (the b-purity shown in Fig. 3a). Other variations change the purities 
very little. 

The differences between the standard results and those found using each of the above condi- 
tions are used to define symmetric systematic uncertainties. To reduce the influence of statistical 
fluctuations, the systematic uncertainties from all sources are determined for a few larger Qj Ct 
bins, each of them exactly covering two or three original bins. The systematic deviation found 
for this larger bin is then assigned to all original bins contained in it. The total systematic 
uncertainty is defined as the quadratic sum of these deviations. 

6 Monte Carlo comparison of biased and unbiased jets 

As discussed above, jets found using a jet algorithm are biased and in this sense are less suit- 
able for comparison with theory than unbiased jets. To assess the difference between biased 
and unbiased jets, we perform a comparison of their properties using hadron level MC event 
samples. For this purpose, we choose HERWIG because it contains an event generator for gg 
events from a colour singlet point source and because it describes well gluon jet properties [4]. 

The conclusions from the comparison of the biased and unbiased jet fragmentation functions 
are basically independent of the Monte Carlo model and jet algorithm used in the analysis, 
therefore, as an example, we show in Fig. 6 the comparison for HERWIG 6.2 and the Durham 
jet algorithm. The results correspond to the hadron level described in Section 5.4. The three-jet 
events (i.e. containing biased jets) are generated at y/s = 91.2 GeV. The inclusive hadronic 
events (no jet finder used, so containing unbiased jets) are generated separately at values of 
yt~s corresponding to twice the central values of Q )et in the individual Q )et intervals used in 
the analysis of three-jet events. Differences between biased and unbiased jet properties are 
expected due to different scales used (<5jet vs - and different number of jets per event 

(two hemispheres vs. three jets found by a jet algorithm and spatially restricted by the minimum 
inter-jet angle of 30°). We point out four regions of phase space where the differences between 
the biased and unbiased jet fragmentation functions are larger than 15%: 

a) Small scales with small xe for all fragmentation functions: This difference, which decreases 

with increasing scale and xe, may in part be explained by hadron mass effect. At small 
c.m.s. energies, hadron masses are not negligible with respect to jet energies, causing a 
suppression of the fragmentation functions at very low xe- This effect is not present in 
theory (hadrons are taken to be massless) and is less strong in three-jet events (the mean 
value of .Ejet in the first Qj Ct bin is about 13 GeV) so one can expect the three-jet data 
to be better described by theory than the unbiased jet data in these bins. 

b) Small scales with large xe for b jet fragmentation functions: Since this difference increases 

with increasing xe and decreasing scale, it might be explained by the b-quark mass effect, 
i.e. by the ratio m\>/ E^ et . At small c.m.s. energies, just above the bb production threshold 
(y/s = 2m h pa 10 GeV), the above ratio is close to 100% and almost all particles picked up 
in the hemispheres come from decays of B hadrons, resulting in a very small probability 
to produce a particle with xe close to unity. As the scale increases, there is more and 
more phase space for gluon showers, leading to more hadrons from the fragmentation 
process. However, the number of radiated gluons is limited by the so-called "dead cone 
effect" [46], i.e. by a suppression of the gluon emission within an angle of order m-^/ E- ]et . 
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In three-jet events, the ratio m^/E- iet starts at a much smaller value than in hemisphere 
events (since the mean jet energy in the first Qj e t bin is about 13 GeV) leading to much 
more phase space for gluon showering compared to hemisphere jets with the same value 
of scale. In QCD calculations based on unbiased jets, this ratio can be identified with 
mass terms of the type m q /Q where Q is some hard scale. In current NLO calculations, 
these mass terms are not considered. As will be seen later, the three-jet data and theory 
behave similarly in the region of small scales. This similarity suggests that missing mass 
terms in theory may behave like m^/E^. 

c) Large xe for gluon jet fragmentation functions: The sizable discrepancy observed for x E > 

0.6 clearly suggests a bias in the gluon jet results. It appears to be more appropriate [47] 
to consider for example both the energy scale and the exact virtuality scale and to boost 
to a frame in which the two scales are equal. MC studies recently presented by OPAL 
in [14] demonstrate that such boosted gluon jets are less biased than those from our study, 
in particular in the regions of very small and large Xe- 

d) The last scale bin for all quark jet fragmentation functions: The observed difference is 

larger than 15% in the xe ranges of 0.01-0.07 and 0.40-0.90. Although biased jets in the 
interval 44 < Q- ]Ct < 46 GeV should in principle resemble hemispheres of the same energy 
(due to large angles $ reaching up to 165°), we found that the soft particle multiplicity 
differs between the two cases. Therefore this difference is considered to represent a true 
bias of biased jets. 

The comparisons made in this MC study suggest that biased jets are less sensitive to hadron 
and b-quark mass effects than unbiased jets. This implies that biased jets tend to be more 
appropriate for comparisons with theory than unbiased jets in the regions of low scale with low 
xe, and in case of b jets, also at low scale with high xe- 

7 NLO predictions 

The results are compared to theoretical predictions by three groups, namely Kniehl, Kramer and 
Potter (KKP) [17], Kretzer (Kr) [18] and Bourhis, Fontannaz, Guillet and Werlen (BFGW) [19]. 
The three groups provide numerical values of the quantity defined in Eq. (2), up to the next- 
to- leading order in as- This means that in the extraction of these predictions from measured 
charged particle momentum distributions, the hard scattering cross section for the production 
of a parton in e + e~ annihilation is evaluated to an accuracy of the order as, while the splitting 
functions describing the scale dependence are evaluated to an accuracy of the order a 2 s . We 
stress that these NLO predictions correspond to an unbiased jet definition. The scale evolu- 
tion via DGLAP evolution equations is performed starting from fragmentation functions at 
a fixed input scale, extracted from existing measurements. In each of these calculations, the 
renormalization and fragmentation scales are set equal to the hard scale Q. The calculations, 
nevertheless, differ in a number of important aspects, such as the choice of data sets, the def- 
inition of the scale Q, the fit ranges, the prescription for the number of active flavours in the 
evolution of fragmentation functions and partonic cross sections, and the treatment of heavy 
flavours and gluons. 

More specifically, in [17] the evolution of the b jet fragmentation function starts at scale 
Q = 2m b where m b is the b-quark mass put equal to 4.5 GeV. The number of active flavours, 
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^active; i s driven by twice the quark mass, 2m q (^ ctive = 4 for 2m c < Q < 2m b and similarly for 
other flavours). The QCD scale parameter for five flavours and the MS renormalization scheme, 
A^g-, is set equal to 0.213 GeV. In [18] the start of the b jet fragmentation function evolution is 
at the scale Q = m b , iV/ ctive is driven by 2m q and = 0.168 GeV. In [19] the fragmentation 
functions are evolved using an "optimal" scale, Q op t, given by the relation Q -jgf |Q=Q opt = 0. 
The evolution of the b jet fragmentation function starts at scale Q = m b , and A^ =0.300 GeV. 

The predictions for quark jet fragmentation functions by KKP, Kr and BFGW were made 
using data from [8-10,48] or similar results. Concerning the predictions for gluon jet fragmen- 
tation functions, it is important to note that in [17] a fit was made to the unbiased [4] and 
biased [49] jet data, in [18] the predictions were obtained from the evolution and the NLO 
correction to the e + e~ cross section and in [19] a fit was made to large p T charged particle 
data [50]. Therefore, the experimental input for gluon jets is very different in the three calcu- 
lations. The fit ranges used by KKP, Kr and BFGW were 0.1 < xe < 1.0, 0.05 < xe < 0.8 and 
0.12 < xe < 0.9, respectively. We obtained the NLO predictions of Kr and BFGW using the 
code [51] where they are provided in parameterised forms. The relative difference between the 
parameterisation and the exact evolution for predictions by Kr are smaller than 3% and 10% 
for xe < 0.75 and x E < 0.90, respectively. All the NLO curves by KKP shown in this analysis 
correspond to the exact scale evolution. 

We point out that in the NLO predictions, the NLO (of the order as) corrections to the hard 
subprocess correspond to inclusive hadron production. For three-jet events, NLO corrections 
are not available and are expected to depend on the jet algorithm used. Our assumption in 
this analysis is that where the biased jet data are observed to be in a good agreement with the 
unbiased jet data, the unknown NLO corrections are apparently small, and the biased jet results 
can be compared to the existing NLO predictions. Despite the sizable differences between the 
biased and unbiased jet MC results reported in points a) and b) of Section 6, the biased jet 
data at low scales are still considered to be appropriate for such a comparison for the reasons 
mentioned at the end of Section 6. 

8 Results 

In the following, the results from this analysis are compared with existing measurements as 
well as with various fragmentation models and theoretical NLO predictions. The fragmentation 
functions are presented either with emphasis on the scale dependence or the xe dependence. 
The scale dependent fragmentation functions are plotted in several xe intervals as functions of 
scale. For a given bin of scale, the data or MC point is placed at the value of the scale at which 
the NLO prediction is equal to its mean value over this bin [52]. An analogous prescription 
is applied for the xe dependent fragmentation functions. Since in the following, the biased 
and unbiased jet results are often plotted on the same figure, we have to accommodate the 
differences between scale definitions and number of jets from which the fragmentation functions 
were extracted. Therefore the term scale in the following figures stands for Qj et in case of biased 
jets and \fsj2 in case of unbiased jets. The published unbiased jet results are scaled by \ since 
they refer to the entire event, thus to two jets. For the NLO predictions, the same prescription 
as for the published unbiased jet data is applied. 
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8.1 Scale dependence 

In Figs. 7-10 and in Tables 3-6 the results for the udsc, b, gluon and flavour inclusive jet frag- 
mentation functions are presented. The LEP1 unbiased jet data correspond to y/s = 91.2 GeV. 
Concerning the LEP2 unbiased jets, the b jet fragmentation functions are measured in the entire 
available y/s range of 183-209 GeV. The corresponding data points are placed at (y/s) = 197 
GeV, where (y/s) is the luminosity weighted value of y/s. The udsc and flavour inclusive jet frag- 
mentation functions are measured in three y/s intervals: 183-189, 192-202 and 204-209 GeV. 
The corresponding data points are placed at (y/s) = 187.6, 198.0 and 206.2 GeV, respectively. 
The quark biased jet data from LEP1 cover the region Qj Ct = 4-42 GeV, while those from 
LEP2 cover the region Q )et = 30-105 GeV. The results from the region 0.01 < x E < 0.03 are 
not shown but they are discussed in Section 8.2. The results are found to be consistent with 
previous measurements. The fragmentation functions from unbiased quark jets agree to within 
the total uncertainties with previous OPAL unbiased jet measurements of flavour inclusive and 
b jet fragmentation functions at y/s = 91.2 GeV in [10] and flavour inclusive jet fragmentation 
functions at y/s = 192-209 GeV in [35] (not shown). Similarly, the udsc and gluon fragmen- 
tation functions from biased jets agree with similar measurements presented by the DELPHI 
Collaboration [5] for Qj Ct scales between 4 and 30 GeV (not shown). Finally, our gluon jet 
results are seen to be consistent with the results of the g inc i jets [4] at 40.1 GeV, see Fig. 9. 
The other results from our study represent first measurements, specifically the udsc jet results 
above 45.6 GeV, the gluon jet results above 30 GeV (except for the gi nc i jets), and the b jet 
results at all scales except 45.6 GeV. 

The data are compared to the theoretical predictions described in Section 7. For the udsc 
jet fragmentation function (Fig. 7), all three theoretical predictions give a good description 
in the entire measured phase space, except for the lowest x E bin where the KKP calculations 
overestimate the data, and the highest x E bin where the data are underestimated by the Kr 
and BFGW calculations. 

The situation is rather different for the b and gluon jet fragmentation functions (Figs. 8 
and 9) where the description of the data by the NLO predictions is worse and where there are 
significant differences between individual NLO results. The latter is, nevertheless, expected 
due to differences in the calculations such as those discussed in Section 7. In Fig. 8 the KKP 
prediction is deficient with respect to the data for x E > 0.12. As shown in [9], with rising 
particle momentum, this region is increasingly populated by the products of B hadron decays. 
It is, however, important to note that these B hadron decay products are indirectly included 
in theory predictions since they are present in the data sets to which the fits were made. 

For the gluon jet fragmentation functions, the two alternative methods of identifying gluon 
jets described in Section 5.2 are examined, see Fig. 9 and Table 5. The Q )Ct binning is not 
the same for the two methods because of their different regions of applicability. In the LEP1 
samples, the interval Q )et = 4-42 GeV is used for the b-tag method, while for the energy- 
ordering method, the Qj et spectra of jets 2 and 3 overlap in the interval Qj et = 6-27 GeV as 
mentioned in Section 5.2.2. In the LEP2 samples, the results correspond to the interval Qj et = 
30-70 GeV for the b-tag method, where only jets 2 and 3 are used, and to the interval Qj et = 
30-60 GeV for the energy-ordering method. A satisfactory correspondence between the b-tag 
and energy-ordering methods is found in the entire scale range accessible. The data tend to 
show larger scaling violations than predicted by any of the calculations. 
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The results for the flavour inclusive jet fragmentation functions are presented in Fig. 10 and 
in Table 6. The results are compared with published unbiased jet data from lower energy e + e~ 
experiments (TASSO, MARK II, TPC and AMY) [8] and previous OPAL results [32-34]. We 
note that the fragmentation functions measured by TASSO, MARK II and AMY are defined 
via x p = 2p/ y/s, where p is particle momentum, rather than via x E used in the present analysis. 
This difference in definition leads to non- negligible differences in the region of xe < 0.1 and 
y/s < 22 GeV, therefore the published data from this region are not shown in Fig. 10. The 
results from the current study are seen to be consistent with the previous results. The data 
are also compared to the NLO predictions of KKP, Kr and BFGW. All three predictions give a 
reasonable description of the data in the central region of x E (0.06 < x E < 0.60) and over the 
entire scale range. 

A good correspondence is found between the results from biased and unbiased jets in all 
four figures. This observation suggests that Qj et is an appropriate choice of scale in three-jet 
events with a general topology. A similar conclusion was previously presented in [5]. The Monte 
Carlo study described in Section 6, however, demonstrates that the bias introduced by using 
jet algorithms in the gluon jet identification is not negligible for xe > 0.6. In each of these 
figures, the scaling violation seen in the data is positive for low xe and negative for high xe- It 
is more pronounced in the gluon jets than in the quark jets. 



8.2 a?E-dependence 

In Section 6, we noted the region of small xe with small scales where large differences between 
biased and unbiased jet fragmentation functions constructed from hadron level MC were ob- 
served. In Fig. 11, this observation is confronted with data. We plot again the unbiased jet 
data of TASSO and the biased jet data from our analysis (Table 6), the latter in those Qj Ct 
bins which correspond well to the c.m.s. energies used in TASSO measurement. We transforme 
x p — > xe using the pion mass and shifte the TASSO points accordingly. In the first scale 
bin, the unbiased jet fragmentation functions exhibit turn-over points at very low xe, while 
the biased jet data grow steeply with decreasing xe- This difference qualitatively confirms the 
observation we made in point a) of Section 6 using MC jet samples. Further, as anticipated in 
Section 6, the biased jet data agree better with theory than the unbiased jet data. 

In Figs. 12-17 we present the results shown in Figs. 7-9 but now in a finer x E binning and 
with the additional data from the region 0.01 < xe < 0.03, see Table 7. We integrate over 
the fragmentation functions in four or five scale intervals: Qj Ct = 4-9, 9-19, 19-30, 30-70 and 
y/s/ 2 = 91.5-104.5 GeV. Reference values for these intervals, evaluated as explained at the 
beginning of this section, are 6.4, 13.4, 24.0, 46.5 (48.5 for gluons) and 98.5 GeV, respectively. 
In the lowest scale interval, the data in the region of 0.01 < xe < 0.02 are not measured due 
to the large dependence on the jet algorithm. 

In Figs. 12-14 the data are compared to the NLO predictions. In general, the theory pre- 
dictions are in a good agreement with the measurements of the udsc jet fragmentation function 
(Fig. 12). We observe that the data in the region of low x E are overestimated by the predictions 
of KKP, while they are in agreement with those of Kr and BFGW. For high x E , the data prefer 
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the KKP predictions but the differences between the predictions decrease with increasing scale. 
In Fig. 13 the measured b jet fragmentation function is shown together with the published re- 
sults from DELPHI [9] and TPC [48]. Analogously to Fig. 8, the spread of the NLO predictions 
is larger than that for the udsc jet fragmentation functions. The NLO predictions by Kr are 
seen to provide a reasonable description of all the b-jet data, while those by KKP and BFGW 
generally overestimate the data in the region of low x E and underestimate them for large xe- 
A possible explanation for this difference is that, unlike KKP or BFGW, the fitting procedure 
of Kr includes both the low x^ (down to x^ = 0.05) and low scale data (TPC data [48] taken 
at t/s =29 GeV). In Fig. 14 the measured gluon jet fragmentation functions are shown along 
with the OPAL [14] measurement at E- )Ct = 14.24 GeV. An overall agreement is found between 
the results of the boost method and the method used here. The observed sizable spread of the 
NLO predictions is expected because of the different approaches to the fitting procedures of 
the gluon jet data (see Section 7). 

To test various fragmentation models, the data are also compared in Figs. 15-17 to the 
hadron level predictions of the PYTHIA 6.125, HERWIG 6.2 and ARIADNE 4.08 MC event 
generators. The hadron level is defined in Section 5.4. Globally, all MC models give a more 
satisfactory description of the data than do the NLO predictions. This is presumably due to the 
fact that the biased jet data are compared to the biased jet MC predictions and the unbiased 
jet data to the unbiased jet MC predictions. We note that although all MC models used in this 
study were previously tuned to LEP1 data, they still provide a good description of the LEP2 
data. There exist some discrepancies in the description of the gluon jet data in the region of 
high xe with small scales (Fig. 17). A good agreement is achieved for the b jet fragmentation 
functions by all three models. 



8.3 Charged particle multiplicities 



By integrating the unbiased jet fragmentation functions, the charged particle multiplicities in 
udsc, b and inclusive hadronic events can be obtained. The results for the LEP2 data are 
presented in the y/s intervals specified above, namely 183-189, 192-202 and 204-209 GeV for 
the inclusive hadronic and udsc events and 183-209 GeV for the b events. 



(V5> 
91.2 GeV 

187.6 GeV 

197.0 GeV 

198.0 GeV 

206.2 GeV 



in 



incl\ 
ch / 



20.93 ±0.01 ±0.23 
26.80 ±0.24 ±0.46 

27.68 ±0.26 ±0.50 
27.75 ±0.29 ±0.67 



«h dsc ) 
20.32 ±0.03 ±0.27 

26.43 ±0.26 ±0.81 

27.38 ±0.31 ±0.85 
26.87 ±0.32 ±0.99 



23.28 ±0.09 ±0.70 
30.01 ±0.53 ±0.82 



The results are found to be in agreement with the previous unbiased jet measurements [9, 10, 
53]. We also observed a good agreement between the data and predictions of the three MC 
models used in this analysis. The results for (n^ sc ) above 91.2 GeV energy represent new 
measurements. 



9 Conclusions 

Scaling violations of quark and gluon jet fragmentation functions are studied in e + e~ annihila- 
tions at \/s = 91.2 and 183-209 GeV using data collected with the OPAL detector at LEP. The 



19 



scale dependence of the flavour inclusive, udsc and b fragmentation functions from unbiased 
jets is measured at \fs/2 = 45.6 and 91.5-104.5 GeV. Biased jets are used to extract the flavour 
inclusive, udsc and b, and gluon fragmentation functions in the ranges Qj e t — 4-42, 4-105 and 
4-70 GeV, respectively, where Qj e t is the jet energy scale. Three methods are used to extract 
the fragmentation functions, namely the b-tag and energy-ordering methods for biased jets, 
and the hemisphere method for unbiased jets. The results obtained using these methods are 
found to be consistent with each other. The udsc jet results above the scale of 45.6 GeV, the 
gluon jet results above 30 GeV (except for the scale of 40.1 GeV), and the b jet results at all 
scales except 45.6 GeV represent new measurements. The results of this analysis are compared 
with existing lower energy e + e~ data and with previous results from DELPHI and OPAL. The 
overall consistency of the biased jet results with the unbiased jet results suggests that Qjet is 
a generally appropriate scale in events with a general three-jet topology. The scaling violation 
is observed to be positive for lower xe and negative for higher xe, for all the types of fragmen- 
tation functions. The gluon jet fragmentation function exhibits stronger scaling violation than 
that of udsc jets. 

The bias of the procedure used to construct biased jet fragmentation functions is estimated 
by studying hadron level Monte Carlo generator events. In explaining the observed differences 
between biased and unbiased jet results, we note the effects of non-negligible masses of hadrons 
and b-quarks at low scales. Due to the considerable bias found for the gluon jet fragmentation 
functions in the region of xe > 0.6, precautions should be taken when comparing the biased 
gluon jet results with theory. 

The data are compared to the predictions of NLO calculations. In a wide range of xe, 
all calculations satisfactorily describe the data for the udsc jet fragmentation functions. The 
description is worse and the spread between the predictions larger for the b and gluon jet frag- 
mentation functions, in particular in regions of very low and high x E . 

The data are also compared with predictions of three Monte Carlo models, PYTHIA 6.125, 
HERWIG 6.2 and ARIADNE 4.08. A reasonable agreement with data is observed for all mod- 
els, except for high xe region with small scales (< 14 GeV) in case of the udsc and gluon jet 
fragmentation functions. 

The charged particle multiplicities of udsc, b and inclusive hadronic events are obtained 
by integrating the measured fragmentation functions. All values are found to be in agreement 
with previous measurements, where available. 
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Selection 


Data LEP1 


Data LEP2 


BG(LEP2) 


Hadronic events 


2 387 227 


10 866 (12 653) 


11% (14%) 


three-jet events 


965 513 


6 177 


16% 


udsc jets 


2 675 679 


16 344 


16% 


b-tag jets 


83 549 


820 


9% 


Gluon jets 


73 620 


729 


9% 


udsc hemispheres 


4 740 774 


20 146 


11% 


b-tag hemispheres 


33 680 


1 586 


5% 



Table 1: Statistics of the LEP1 and LEP2 data samples. The number of hadronic events is 
given (the numbers in brackets correspond to LEP2 events used for three-jet analysis) and the 
characteristics of the biased jet samples (number of three-jet events, udsc, b-tag and gluon jets) 
and the unbiased jet samples (number of udsc and b-tag hemispheres) are shown. The jets are 
found by the Durham jet algorithm. Also indicated is the percentage of the remaining 4-fermion 
background (BG) for the LEP2 data. For the LEP1 data, the background is negligible. 



Cuts 


Loss 


[%] 




LEP1 


LEP2 


Particle multiplicity per jet > 2 


0.7 


1.7 


Sum of inter-jet angles > 358° 


3.9 


2.3 


Polar jet angle cos0 jet | < 0.90 (0.95) 


8.4 


2.3 


Corrected jet energy > 5 GeV 


11.2 


5.9 


Inter-jet angle > 30° 


43.3 


43.2 



Table 2: Jet selection cuts for the LEP1 and LEP2 data and the reduction of statistics found 
by imposing each cut individually after the hadronic event selection. The jets are found by the 
Durham jet algorithm. The cut value given in brackets corresponds to the LEP2 selection. 
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Table 3: The udsc jet fragmentation function in bins of xe and scale. The scale denotes Qjct 
for the biased jets and is given by the intervals, while it denotes y/s/2 for the unbiased jets and 
is given by the single values. The first uncertainty is statistical, the second systematic. These 
data are displayed in Fig. 7. 
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Table 4: The b jet fragmentation function in bins of xe and scale. The scale denotes Qj et for the 
biased jets and is given by the intervals, while it denotes \fs/2 for the unbiased jets and is given 
by the single values. The first uncertainty is statistical, the second systematic. These data are 
displayed in Fig. 8. In the region 0.48 < x E < 0.90 and Qjet = 30-70 GeV, no measurement 
was possible due to low statistics. 
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1 djV ch 



(EO) 
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Table 5: The gluon jet fragmentation functions in bins of x E and scale Qjct obtained from 
the biased jets using the b-tag method (BT) and the energy-ordering method (EO). The first 
uncertainty is statistical, the second systematic. These data are displayed in Fig. 9. In the 
region 0.48 < x E < 0.90 and Q jct = 30-42 GeV for the b-tag method and Q jct = 24-60 GeV 
for the energy-ordering method, no measurement was possible due to low statistics. 
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53.38 
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±0.11 

±0.055 

±0.046 
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12.0 
15.0 
19.0 
24.0 
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- 6.5 

- 9.0 

- 12.0 

- 15.0 

- 19.0 

- 24.0 

- 30.0 

- 42.0 
45.6 
93.8 
99.0 
103.1 



30.0 

31.6 

28.2 

29.3 

30.2 

29.79 

30.11 

30.26 

32.32 

31.56 

31.93 

31.70 



±2.8 

±1.1 

±0.5 

±0.5 

±0.4 

±0.35 

±0.32 

±0.15 

±0.03 

±0.38 

±0.43 

±0.46 



±2.8 
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±1.0 
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0.664 

0.658 

0.625 

0.575 

0.564 

0.546 

0.507 

0.508 
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Table 6: The flavour inclusive jet fragmentation functions in bins of xe and scale. The scale denotes 
Qjct for the biased jets and is given by the intervals, while it denotes y/s/2 for the unbiased jets and 
is given by the single values. The first uncertainty is statistical, the second systematic. These data 
are displayed in Figs. 10 and 11. 
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Table 7: The udsc-quark, b-quark and gluon fragmentation functions in bins of xe and scale. The 
scale denotes Qj e t for the biased jets (the first four intervals) and \/s/2 for the unbiased jets (the single 
value and the last interval). The first uncertainty is statistical, the second systematic. These data are 
displayed in Figs. 12-17. 
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Figure 1: (a) Jet energy and (b) jet scale Qjet distributions for the selected three-jet events in the 
LEP1 sample. The solid histograms represent the JETSET 7.4 and the dotted histograms the 
HERWIG 6.2 predictions using the Durham jet algorithm. The data are shown with statistical 
uncertainties only. 
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Figure 2: (a) Jet energy and (b) jet scale Qj et distributions for the selected three-jet events in 
the LEP2 sample. The solid histograms represent the sum of the PYTHIA 6.125 and back- 
ground (BG) predictions, the dotted histograms the sum of the HERWIG 6.2 and background 
predictions and the shaded histograms the prediction of the model GRC4F for the background, 
all using the Durham jet algorithm. The data are shown with statistical uncertainties only. 



31 



1 

0.8 
0.6 
0.4 
0.2 




(a) b-tag jets (y/s = 91.2 GeV) 



O O O O O O o 

• Purity 
o Efficiency 



o o 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

VNN 




Figure 3: The purity (solid circles) and efficiency (open circles) as a function of the neural 
network output VNN for (a) the LEP1 b-tag jet sample and (b) the LEP2 b-tag jet sample, 
obtained from the JETSET 7.4 for the LEP1 events and PYTHIA 6.125 for the LEP2 events 
using the Durham jet algorithm. 
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Figure 4: The purity matrix as a function of Qj e t scale for the LEP1 (open circles) and the 
LEP2 (solid circles) three-jet events, obtained from the JETSET 7.4 for the LEP1 events and 
PYTHIA 6.125 for the LEP2 events using the Durham jet algorithm. The flavour assignment 
is used for matching the primary outgoing partons to hadron jets (see text). 
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Figure 5: The quark purities of jet 1 and 2 samples (Pi q and P 2q ) and the gluon purity of jet 3 
sample (P3 g ) as a function of the scale Q- ]Ct using the Durham jet algorithm. Either the matrix 
element (ME) information (data as symbols and JETSET 7.4 as the solid histograms) or the 
matching method (JETSET 7.4 as dashed histograms) is used. Only statistical uncertainties 
are shown. 
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Figure 6: The fragmentation functions in bins of Xe and scale as obtained from hadron level 
events generated with HERWIG 6.2. The scale stands for Qj et in the case of biased jets (gen- 
erated at y/s = 91.2 GeV) and for \fsj1 in the case of unbiased jets. The Durham algorithm 
is used to find jets. The statistical uncertainties are smaller than the symbols. 
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Figure 7: Scale dependence of the udsc jet fragmentation functions in different bins. The 
scale denotes Qj Ct for the biased jets and \fsj1 for the unbiased jets. The inner error bars 
indicate the statistical uncertainties, the total error bars show the statistical and systematic 
uncertainties added in quadrature. The values are given in Table 3. The data are compared to 
the NLO predictions by KKP [17], Kr [18] and BFGW [19]. 
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Figure 8: Scale dependence of the b jet fragmentation functions in different xe bins. The scale 
denotes Q- )et for the biased jets and \fsj2 for the unbiased jets. The inner error bars indicate the 
statistical uncertainties, the total error bars show the statistical and systematic uncertainties 
added in quadrature. The values are given in Table 4. The data are compared to the NLO 
predictions by KKP [17], Kr [18] and BFGW [19]. 
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Figure 9: Scale dependence of the gluon jet fragmentation functions in different xe bins. The 
scale denotes Q- )Ct for the biased jets and E- )Ct for the published g inc i jets (OPAL [4]). The 
results from the biased jets using the b-tag (BT) and the energy-ordering method (EO) are 
shown. The inner error bars indicate the statistical uncertainties, the total error bars show the 
statistical and systematic uncertainties added in quadrature. The values are given in Table 5. 
The data are compared to the NLO predictions by KKP [17], Kr [18] and BFGW [19]. 
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Figure 10: Scale dependence of the flavour inclusive fragmentation functions in different xe 
bins. The scale denotes Q ]et for the biased jets and y/s/2 for the unbiased jets. The inner 
error bars indicate the statistical uncertainties, the total error bars show the statistical and 
systematic uncertainties added in quadrature. The values are given in Table 6. In addition, 
the published unbiased jet data by TASSO, TPC, MARK II, AMY [8] and OPAL [32-34] are 
shown. The data are compared to the NLO predictions by KKP [17], Kr [18] and BFGW [19]. 
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Figure 11: x-^ dependence of the flavour inclusive fragmentation functions at different scales. 
The inner error bars indicate the statistical uncertainties, the total error bars show the statistical 
and systematic uncertainties added in quadrature. The values are given in Table 6. The results 
from the biased jets are compared to the results from the unbiased jets of TASSO [8]. The data 
are compared to the NLO predictions by KKP [17], Kr [18] and BFGW [19]. The results are 
scaled as indicated in figure. 
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Figure 12: xe dependence of the udsc jet fragmentation functions at different scales. The inner 
error bars indicate the statistical uncertainties, the total error bars show the statistical and 
systematic uncertainties added in quadrature. The values are given in Table 7. The data are 
compared to the NLO predictions by KKP [17], Kr [18] and BFGW [19]. The results are scaled 
as indicated in figure. 
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Figure 13: dependence of the b jet fragmentation functions at different scales. The inner 
error bars indicate the statistical uncertainties, the total error bars show the statistical and 
systematic uncertainties added in quadrature. The values are given in Table 7. For comparison, 
the published results on the unbiased jets of DELPHI [9] and the results based on TPC data [8] 
are shown. The data are compared to the NLO predictions by KKP [17], Kr [18] and BFGW 
[19]. The results are scaled as indicated in figure. 
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Figure 14: xe dependence of the gluon jet fragmentation functions at different scales. The 
inner error bars indicate the statistical uncertainties, the total error bars show the statistical 
and systematic uncertainties added in quadrature. The values are given in Table 7. Also shown 
are the recent OPAL results on the boosted gluon jets [14]. The data are compared to the NLO 
predictions by KKP [17], Kr [18] and BFGW [19]. The results are scaled as indicated in figure. 
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Figure 15: xe dependence of the udsc jet fragmentation functions at different scales compared 
with the PYTHIA 6.125, HERWIG 6.2 and ARIADNE 4.08 Monte Carlo predictions. The 
inner error bars indicate the statistical uncertainties, the total error bars show the statistical 
and systematic uncertainties added in quadrature. The values are given in Table 7. The results 
are scaled as indicated in figure. 
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Figure 16: xe dependence of the b jet fragmentation functions at different scales compared 
with the PYTHIA 6.125, HERWIG 6.2 and ARIADNE 4.08 Monte Carlo predictions. The 
inner error bars indicate the statistical uncertainties, the total error bars show the statistical 
and systematic uncertainties added in quadrature. The values are given in Table 7. For 
comparison, the published results on the unbiased jets of DELPHI [9] are shown. The results 
are scaled as indicated in figure. 
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Figure 17: x-^ dependence of the gluon jet fragmentation functions at different scales compared 
with the PYTHIA 6.125, HERWIG 6.2 and ARIADNE 4.08 Monte Carlo predictions. The 
inner error bars indicate the statistical uncertainties, the total error bars show the statistical 
and systematic uncertainties added in quadrature. The values are given in Table 7. The results 
are scaled as indicated in figure. 
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